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Abstract 

An influence of a strong external magnetic field on the neutrino self-energy 
operator is investigated. The width of the neutrino decay into the electron 
and W boson, and the mean free path of an ultra-high energy neutrino in a 
strong magnetic field are calculated. A kind of energy cutoff for neutrinos 
propagating in a strong field is defined. 
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1. Introduction 



The solution to the solar-neutrino puzzle in a unique experiment at the 
heavy-water detector installed at the Sudbury Neutrino Observatory has un- 
doubtedly been the most important achievement of neutrino physics within 
the last decades. This experiment confirmed B. Pontecorvo's key idea con- 
cerning neutrino oscillations and, along with experiments that studied at- 
mospheric and reactor neutrinos, thereby proved the existence of a nonzero 
neutrino mass and the existence of mixing in the lepton sector, see e.g. [l[ 
and the references cited therein. In this connection, the problem of study- 
ing the possible effect of an active environment, including a strong magnetic 
field, on the dispersion properties of the neutrinos becomes quite important. 

An analysis of the effect of an external medium on neutrino properties 
relies on calculating the neutrino self-energy operator E(p), from which one 
can extract the neutrino dispersion relation, and in part the imaginary part 
of the neutrino self-energy in medium, defining the width of the neutrino 
decay into the W + boson and a charged lepton, v — > £~W + . In this paper, 
we consider an electron as a charged lepton, but all the formulas are valid to 
the muon and r lepton as well. 

A literature search reveals that calculations of the neutrino dispersion 
relation in external magnetic fields have a long history 0, 0, 0, 0, @, Q, 0, @]. 
To compare the different results we analyse the neutrino self-energy operator 
that is defined in terms of the invariant amplitude for the transition 
v e — f v e by the relation 

M[v e u e ) = - Hp) E(p) u(p)] = -Tr [£(p) p{p)} , (1) 

where p = (E, p) is the neutrino four-momentum, p(p) = v{p)v{p) is the 
neutrino density matrix. On the other hand, the additional energy AE ac- 
quired by a neutrino in an external magnetic field is defined via the invariant 
amplitude ([TJ as follows: 

*E = ~M{u e ^v e ). (2) 

The S matrix element for the transition v e — > v e corresponds to the Feyn- 
man diagrams shown in Fig. [1] where double lines denote exact propagators 
in the presence of an external magnetic field. A detailed description of the 
calculational techniques for the neutrino self-energy operator in exter- 
nal electromagnetic fields can be found e.g. in Ref. see also 0, 0, liol ]. 
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The relevant iS-matrix element can be used to deduce, in a standard way, 
the invariant amplitude flTJ, whence the neutrino self-energy operator takes 
the form 



+ — — (m e R — m v L) J^\p) (m e L — m v R) . (3) 

Here, g is the Standard Model electroweak coupling constant; 7 Q are the 
Dirac matrices; L = (1 — 7s)/2 and i? = (1 + 7 5 )/2 are, respectively, the left- 
and the right-hand projection operator. The integrals introduced in (J3]) have 
the form 

■Cm = / (gj^e^ta-p), 

J< * >(p) = j WT lS{ ' ,)Dm(q ' p) ' (4) 

where S'(g), (q — p) and D^\q — p) are the Fourier transforms of the 
translation-invariant parts of the propagators for the electron, the W~ boson, 
and the charged scalar $ boson, respectively. We note that the quantity 
m„ in ([3]) is in general the nondiagonal Dirac neutrino mass matrix with 
allowance for mixing in the lepton sector. 

The general Lorentz structure of the operator in a magnetic field, 
defined in Eq. ([3]), can be represented in the form \j\ 

S(p) = [A L (iry)+B L (p7)\\+C L (p07)] L 
+ [A R {pi) + B R (pj) 11 + Cr (pfry)] i? 
+ 77V [/Ci + i/C 2 (7^7)] • (5) 

The Lorentz indices of four- vectors and tensors within parentheses are con- 
tracted consecutively, e.g. (p<pj) = p a ^p a i3l 13 ■ Further, is the dimension- 
less tensor of the electromagnetic field, normalized to the external 73-field, 
whereas (p is its dual, 



fa/3 



aft 



B ' 

= -£ a f3^u^ V - (6) 
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Finally, in the frame where only an external magnetic field B is present, 
we take the spatial 3-axis to be directed along B. Four- vectors with the 
indices _L and || belong to the Euclidean {1, 2}-subspace and the Minkowski 
{0, 3}-subspace, correspondingly For example, p± = (0,pi,p2,0) and p\\ = 
(po, 0, 0,^)3). For any four- vectors P and Q we use the notation 



The coefficients Ar, Br, Cr, and in §5§ stem from the Feynman 
diagram involving the scalar $ boson, while the coefficients Al, Bl, and Cl 
contain the contributions from both diagrams. We note that the coefficients 
Al, Ar, and /Ci in ([5]) contain an ultraviolet divergence which is removed by 
the vacuum renormalization of the neutrino wave function and mass. 

Using Eqs. ([I]), §2§ and the neutrino additional energy AE in an 
external magnetic field can be written in the form: 



where v = p/E is the neutrino velocity vector, s is the unit vector of the 
doubled neutrino spin, b is the unit vector along the magnetic field direction, 
and b 4i £ are its transversal and longitudinal components with respect to the 
neurino momentum, b = b t + b^. 

In the previous papers, the neutrino self-energy operator (J3J) was calcu- 
lated in different regions of values of the physical parameters, however, the 
list of these considered regions appears not to be comprehensive. Namely, 
the investigated limiting cases were the following: 

i) a weak field case [eB <C ml) 0, [fjj]; 

ii) a moderately strong field case (ml eB f^w) 

iii) the situation where the neutrino transverse momentum p± with respect 
to the magnetic field is rather high, for example, p± > m w or p±_ ^> 



(PQ)\\ 
(PQ)± 
(PQ) 



(P00Q) = P o Q o -P 3 Q 3 , 
(P W Q) = P 1 Q 1 + P 2 Q 2 , 
(PQ) ]{ -(PQ) ± . 



(7) 



AE 




(8) 
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raw, while the magnetic field strength is not too high, eB <C m 2 e) which 
corresponds to the crossed- field approximation (slU, 0, 

There is yet another region of values of the physical parameters that 
requires a dedicated analysis. We mean here the case of the high neutrino 
transverse momentum, when the magnetic field strength is also rather high, 
thus, the crossed-field approximation is not valid. 

This region of parameter values is of importance in connection with prob- 
lems of the physics of magnetars, the pulsars with superstrong surface mag- 
netic fields (B s ~ 10 15 G). In particular, the possibility of detecting cosmic 
neutrinos of ultrahigh energy, ~ 1 PeV or even higher, from magnetars is 
widely discussed (see, for example, [l2|, |l3|). It looks reasonable that the 
process of emission of neutrinos having such energies cannot be described ad- 
equately without taking into account their interaction with a strong magnetic 
field of a magnetar. 



2. Charged-lepton, W- and <fr-boson propagators in a magnetic 
field 

The Fourier transforms of the translation-invariant parts of the exact 
propagators in an external magnetic field, entering into expressions (jl]) can 
be presented in the Fock proper-time formalism in the following form. The 
lepton propagator is [14 



ds 



iq cos (3 s 
x { [(g7)|| +W| 



cos (3s 



(7^7) • a 

- snips 



(97) 



cos (3s 



(9) 



where /3 = eB and m e is the electron mass. 

Similarly, the VT-boson propagator can be written as 

ds 



G prJ (q) 



-ifij, 



X 



i cos/3s 
{<P<p) P o - (w) P a cos 2(3s - y? pCT sin 2(3s 
And finally, for the $-boson propagator one obtains 
£><*>(g) = / dso- i!! " 



(10) 



(11) 
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where we have chosen the Feynman gauge for the W and <3> bosons and have 
introduced the notation (J = e, W) 

ni = a (mJ-gf) + ^^. (12) 

3. The neutrino decay w — > e~W + in an external electromagnetic 
field 

The probability of the neutrino decay v — > e~W + in an external electro- 
magnetic field is one of the most interesting results that can be extracted 
from the neutrino self-energy operator. This probability can be expressed in 
terms of the imaginary part of the amplitude (PQ) with the neutrino self-energy 
operator (jSJ). 

For simplicity, hereafter we neglect the neutrino mass m„, taking the 
density matrix of the left-handed neutrino as p(p) = (jyy) L. One obtains: 

w(v -> e~W + ) = — Im M[y e -> v e ) 
E 

-ilmTr[S(j9)(p 7 )^] = -2 4 Im ^- ( 13 ) 

An analysis of the neutrino decay v — > e~W + in an external field is of 
interest only at ultrahigh neutrino energies. 

In all previous papers the neutrino decay width in an external electro- 
magnetic field was calculated in the crossed field approximation, in which 
case the width is expressed in terms of the dynamical field parameter x an d 
the lepton mass parameter A: 

X = <? FF f 12 , A = 4. (14) 

The particular case of a crossed field is in fact more general than it may 
seem at first glance. Really, the situation is possible when the field dynam- 
ical parameter x of the relativistic particle propagating in a relatively weak 
electromagnetic field, F < B e (where F means the electric and/or magnetic 
field strength), could appear rather high. In this case the field in the particle 
rest frame can exceed essentially the critical value and is very close to the 
crossed field. Even in a magnetic field whose strength is much greater than 
the critical value, the result obtained in a crossed field will correctly describe 
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the leading contribution to the probability of a process in a pure magnetic 
field, provided that x ^ B/B e . In the frame where the field is pure magnetic 
one, the dynamical field parameter takes the form: 

X = ^l3— • ( 15 ) 



m 



A general expression for the decay width can be written in this case in 
the form [3] 



w(v ->• e W 



YItxE 

dvv [2(1 + v){2 + v) + A (1 - v)(2 - v)} dAi(u) 



X ' [v{l-v)Y'z du ' (16) 



where 



oo 

3 



Ai(u) = - [ dt cos (tu + —\ (17) 



7T 

b 

is the Airy function with the argument: 

_ v + \(l-v) 

[Xv(l-v))*/f ^ 

The derivative of the Airy function is expressed via the modified Bessel func- 
tion K v (x) 

dAiW - » K„(\M. (19) 



du y/3 



7T 



Taking in Eq.f lTB]) the limit x, A 1, one obtains the result which can be 
written in terms of the only modified dynamical field parameter 

(20) 

The range for the £ parameter appears to be rather large, < £ <C 1/V% 
while 1/VA > 1. 
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The decay width takes the form 



V2G F (eBp ± ] 



2 



v>(y -»- e~W + ) = v /t F(0 , (21) 



37r F 



where 



F (O = -zH; I d. — A- 2 /3 ( ; ) • (22) 



1 + /2(l+a;) 3 / 2 



We remind that these formulas are valid in the approximation £ <C my/ jm e . 
The range being very wide for the electron, ( C 1.6 x 10 5 , is not too wide 
for the r lepton, £ -C 45. 

The integration in Eq. ( 122|) can be performed exactly to give 

= (l + f ) exp . (23) 

The formulas (j2~Tj) - fl23|) should be compared with the results of Refs. 0, 
0, [sj . It should be mentioned that the decay width w defined in Refs. 0, 0] is 
the same, in the natural system of units, than the absorption coefficient a 
and the damping rate of the neutrino 7 8 
coefficient a presented in Eq.(25) of Ref. 
and (123]) . However, the angular dependence in our formulas is quite different: 
instead of the factor p\/E — Esin 2 6 standing in our Eq. ( 12T|) . there is the 
factor p ± = Esm6 in Eq.(25) of Ref. fl. 

On the other hand, one can see that our result ( J2TT) - ( 1231) surely contra- 
dicts the Eq. (58) of Ref. jsj , where an attempt was made of reinvestigation 
of the process v — > e~W + in the crossed field approximation. The difference 
is the most essential at small values of £, where the result of Ref. |s| appears 
to be strongly underestimated. 

In the earlier paper by Borisov et al. [3!] the calculations of the process 
v — > e~W + width were performed in the two limiting cases of the small and 
large values of the parameter x- I n the limit \ 2 ^ ^ their result can be 
presented in the form 



One can see that the absorption 
5| looks very similar to our Eqs. ( )2T|) 



w 



m e e5sin#exp ( -Vs , (24) 

VSn V eBp ± J 
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and can be reproduced from the general formulas (1211) 

On the other hand, in the limit x> 1 (C>1/ VX) the result of Ref . |3( 
can be written as 

w = y ^J^ F mw e B s i n q ; (25) 

V2 7T 

and can be reproduced from our more general formula (|T6|) . 

A problem of the decay v — > e~W + has a physical meaning only in the 
fields of the pulsar type, where the field strength is of order of the critical 
value ~ 10 13 G. The above formulas for the probability except for Eq. (125]) are 
applicable for relatively weak fields only, B <C 10 13 G. Taking into account 
the discovery of magnetars which are the neutron stars with the fields ~ 
10 14 — 10 15 G, it is interesting to calculate the probability of the process v — > 
e~W + in such fields where the crossed-field approximation is inapplicable. 

Thus, we will use the following hierarchy of the physical parameters: 
p\ ^> m 2 ^ 3> eB 3> m 2 . A general expression for the process v — > e~W + 
probability can be obtained by the substitution of Eq. ([3]) into Eq. ({TBI with 
taking account of Eqs. §§§ - (TTTT) . After calculations which are not difficult 
but rather cumbersome, the process width can be presented in the form 

w(u e-W+) = Gf(6 ^ P± m , (26) 



where Q(rj) is the function depending on the one parameter r\ only: 



oc 



\ 1 f dy (tanhy) 1 / 2 (sinhy) 2 — y tanhy 
V J V 1 ^ 2 (sinh|/) 2 (y — tanhy) 3 / 2 



o 



y tanhy 



r)(y — tanhy) 



.(28) 



The formulas ( 1261) . (1281) present our main results and are valid in a wide 
region of the parameter r\ values, < r\ <C m^/^eB). The function $(77) is 
essentially simplified at large and small values of the argument. 

In the limit rj ^> 1 one obtains: 

^»l)-^\A(?7-0.3), (29) 
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and the error is less than 1 % for 77 > 10. 

The formulas (l2"B"j) . ( )2"9"|) reproduce the probability (12"T|) . ( I2"3"j) . where the 
limit £ 3> 1 should be taken, and F(£ 1) ~ 1. 

In the other limit 77 <C 1 one obtains 

$(7? < 1) ~ exp (-^j (l - l -r) + |jr/) (30) 

and the error is less than 1 % for 77 < 0.5. 

The formulas obtained allow to establish an upper limit on the energy 
spectrum of neutrinos propagating in a strong magnetic field. Let us take 
the typical size R of the region with the strong magnetic field as R ~ 10 
km. If the neutrino mean free path A = 1/w is much less than the field size, 
A <C R, all the neutrinos are decaying inside such the field. For A = 1 km 
<C R, we can find the cutoff energies E c for the neutrino spectrum, depending 
on the magnetic field strength, as follows: 

i) for relatively weak field, B ~ 0.1 B e ~ 4 x 10 12 G, the neutrino mean 
free path can be obtained from Eq. (T5] 



> 4.9 m x / 219 \ 

B .! sm9 6XP \B 0A E 15 sinfly ' 1 ' 

where B 01 = B/(0.1B e ), Ei 5 = S/(10 15 eV), and the cutoff energy corre- 
sponding to A = 1 km, at £> .i = 1, = ir/2, is 

E c ~ 0.4 x 10 17 eV; (32) 

ii) for relatively strong field, B ~ 10-B e ~ 4 x 10 14 G, the neutrino mean 
free path can be obtained from Eqs. (|2^|) . (15D1 : 

3.2 cm / 4.0 \ 

A " R 3/2 ■ ~ eX P ( B F 2 -2n > ( 33 ) 

where B\§ = B/(10B e ), and the cutoff energy corresponding to A = 1 km, 
at B w — 1, 9 — 7r/2, is 

E c ~ 0.6 x 10 15 eV. (34) 

The results obtained show an essential influence of the intense magnetic 
field on the process v — > e~W + width. Despite the exponential character of 
suppression of the width in a strong field, Eq. fl29l) . as well as in a weak field, 
Eq. (l24"j) . the decay width in a strong field is greater in orders of magnitude 
than the one in a weak field, for the same neutrino energy. 
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Figure 1: Feynman diagrams representing the magnetic- field-induced contribution to the 
neutrino self-energy operator in the Feynman gauge. Double lines correspond to the exact 
propagators for the charged lepton, the W boson, and the nonphysical scalar charged $ 
boson in an external magnetic field. 
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